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The crossover in energy level statistics of a quasi- 1- 
dimensional disordered wire as a function of its length L is 
used, in order to derive its averaged localization length, with- 
out magnetic field, in a magnetic field and for moderate spin 
orbit scattering strength. An analytical function of the mag- 
netic field for the local level spacing is obtained, and found to 
be in excellent agreement with the magnetic field dependent 
activation energy, recently measured in low-mobility quasi- 
one- dimensional wires This formula can be used to extract 
directly and accurately the localization length from magne- 
toresistance experiments. In general, the local level spacing 
is shown to be proportional to the excitation gap of a vir- 
tual particle, moving on a compact symmetric space. Pacs- 
numbers: 72.15.Rn,73.20.Fz,02.20.Qs 

In disordered wires quantum interference results in lo- 
calization of all states for arbitrary disorder strength, if 
the wire is infinite It has been discovered that the 

localization length depends on the global symmetry of 
the wire Q): L c = (3irhvSDo, where (5 = 1,2,4, corre- 
sponding to no magnetic field, finite magnetic field, and 
strong spin- orbit scattering or magnetic impurities, re- 
spectively. v(E) is the electronic density of states in the 
wire. Do is the classical diffusion constant of the elec- 
trons in the wire, and S its crossection. This result was 
obtained by calculating the spatial decay of the density 
correlation function for wires whose thickness exceeds the 
mean free path I. Independently, it was obtained by cal- 
culating the transmission probability through thin, few 
channel wires ||. Recently, the doubling of the local- 
ization length was observed in sub-micron thin wires of 
doped GaAs by Khavin, Gershenson and Bogdanov, who 
found a continously decreasing activation energy when 
the magnetic field is increased, saturating indeed at one 
half of its field free value However, the crossover 
as function of moderate symmetry breaking fields defied 
any attempt to study it with nonperturbative methods. 
The only results were based on a heuristic approach by 
Bouchaud H, a kind of semiclassical analysis by Imry 
and Lerner [J?J and numerical studies [§,U . Recently, how- 
ever, Kolesnikov and Efetov |l0| succeeded to tackle that 
complex problem and derived the density- density corre- 
lation function with the supersymmetry method in the 
crossover regime. The result does not scale with a single 
localization length. 

Here, the problem of the crossover is addressed from 
a different perspective. The statistics of discrete en- 
ergy levels of a finite coherent, disordered metal parti- 
cle is an efficient way to characterize its properties [||. 



This can be studied by calculating a disorder averaged 
autocorrelation function between two energies at a dis- 
tance u> in the energy level spectrum. It is an oscil- 
latory function whose amplitude decays with a power 
law when the energy levels in the vicinity of the cen- 
tral energy E are extended. A Gaussian decaying func- 
tion is a strong indication that all states are localized. 
The autocorrelation function of spectral determinants 
(ASD) is defined by C(w) = C(lo)/C(0), C(u) = 
(det(E + uj/2 - H)det(E -w/2-H)).E\&& central en- 
ergy. We consider the hamiltonian of disordered electrons 



H= [p-qA] 2 /2m + y(x), 



(1) 



where q is the electron charge, c the velocity of light. 
V(x) is taken to be a Gaussian distributed random func- 
tion (V(x)) = 0, and (V(x)V{x')) = (AS£/2ttt)(5(x - 
x'), which models randomly distributed, uncorrelated im- 
purities in the sample. A is the mean level spacing, 1/r 
the elastic scattering rate. The vector potential is used 
in the gauge A = (— By, 0,0), where x is the coordinate 
along the wire of length L, y the one in the direction 
perpendicular both to the wire and the magnetic field B, 
which is directed perpendicular to the wire. The angular 
brackets denote averaging over impurities. 

The ASD appears as an intermediate step in the Grass- 
man replica trick and as the compact sector of the super- 
symmetric field theory of disordered systems . It is a 
non-self-averaging quantity. It has been shown recently, 
that the ASD can distinguish between an uncorrelated 
spectrum of localized states and a correlated spectrum 
of extended states in the vicinity of E. This way, the 
metal- insulator crossover as a function of the length of 
a disordered wire has been explored with the ASD. A 
localization length £ c in a moderate magnetic field has 
been derived as the crossover length scale and shown to 
coincide with L c {j3 = 2) 11 1. The ASD also yields qual- 
itative information on the location of localized states in 
a quantum- Hall- system |l2] ]. 

The impurity averaged ASD can be written as a par- 
tition function |12|] 



C(w) = Trexp(-LH [Q]), 



(2) 



where H is an effective Hamiltonian of matrices Q on 
a compact manifold, determined by the symmetries of 
the Hamiltonian H of disordered electrons. Thus, the 
problem reduces to the one of finding the spectrum of 
the effective Hamiltonian H. 

There is a finite gap Eq between the ground state en- 
ergy and the energy of the next excited state of H(w = 0). 
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For a long wire, LEq 3> 1, the ASD becomes, according 
to Eq. (||), C(uj) — exp(~const.Luj 2 / Eq)- This is typical 
for a a spectrum of localized states |l2j . In the other limit 
LEq <C 1, all modes of H do contribute to the trace in 
the partition function Eq. (||) with equal weight, yield- 
ing the correlation function of a spectrum of extended 
states JO]]. Thus, the crossover length £ c ~ 1/Eq can be 
identified with an averaged localization length. 

Being a product of two spectral determinants, 4a - 
component Grassman fields are needed to get the func- 
tional integral representation of the ASD. Here, a = 1, 
when the Hamiltonian is independent of the spin of the 
electrons, and each level is doubly spin degenerate. There 
is one pair of Grassman fields for each determinant in the 
ASD and each pair is composed of a Grassman field and 
its time reversed one, as obtained by complex conjuga- 
tion, a = 2 has to be considered, when the Hamiltonian 
does depend on spin, as for the case with moderately 
strong magnetic impurity or spin- orbit scattering. This 
necessitates the use of a vector of a spinor and the cor- 
responding time reversed one. One notes a global invari- 
ance under rotations of these vectors. 

Averaging over impurities according to Eq. (1), one 
gets an interacting theory of Grassman fields with inter- 
action strength 1/r. It can be decoupled with a Gaussian 
integral over 4a x 4a matrices which preserve the global 
rotational invariance. The Grassman fields can now be 
integrated out exactly. Finding the saddle point of the 
integral over matrices and integrating out longitudinal 
(massive) modes, for w <c 1/r and 1/t » A, the ASD 
reduces to a functional integral over transverse modes Q, 
being elements of respective symmetric spaces. For dis- 
ordered wires whose thickness W is larger than the mean 
free path I but smaller than the magnetic length W < Ib, 
the action reduces to the one of the one-dimensional com- 
pact nonlinear sigma model jl3|;[l4| , 



C{io) 



J]dQ(x)exp(-F[Q]), 



(3) 



F[Q) = a-^Lcu J dxTr 
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d.i: 



TrA 3 Q(a;). 



(4) 



where Leu — Lc{(3 = 2) = 2-kTivSDq is the localization 
length in the wire in a moderately strong magnetic field 
J|. The integral over the crossection S of the wire is 
done, giving y 2 . Here, and in the following, A^ are the 
Pauli matrices in the subbasis of the left and the right 
spectral determinant, r, the ones in the subbasis spanned 
by time reversal and Ui the ones in the space spanned by 
the spinor, for i = 1,2, 3. The fluctuations of the ma- 
trices Q are transverse, Q 2 = 1, and restricted by global 
rotational invariance to Q + = Q. For a = 1 in addition 



Q T C — CQ, where C = —it^- This constrains Q to be an 
element of the group defined on the compact symplectic 
symmetric space Sp(2)/( Sp(l) x Sp(l)). For a moderate 
magnetic field, Q is reduced to a 2 x 2- matrix by the bro- 
ken time reversal symmetry. This reduces the space of Q 
to U(2)/(U{1) x U{1)). For a = 2 the matrix C is, due 
to the time reversal of the spinor, substituted by io^ri 
|l4j |. Both magnetic impurities and spin-orbit scattering 
reduce the Q matrix to unity in spin space. Thus, C has 
effectively the form t\. The condition Q T C = CQ leads 
therefore to a new symmetry class, when the spin sym- 
metry is broken but the time reversal symmetry remains 
intact. This is the case for moderately strong spin-orbit 
scattering. Then, Q are 4x4- matrices on the orthogonal 
symmetric space 0(4)/(0(2) xO(2)) (l5j. With magnetic 
impurities both the spin and time reversal symmetry is 
broken, and the Q- matrices are in the unitary symmet- 
ric space Z7(2) / (Z7(l) x U(l)) as for a moderate magnetic 
field and spin degenerate levels. The difference in the 
prefactor a in the action, Eq.(|^), remains. One can ex- 
tend this approach to other compact symmetric spaces 
with physical realizations, see Ref. |l6| for a complete 
classification. For uj/ A < Leu I L, the spatial variation of 
Q can be neglected and one retains the same ASD as for 
random matrices of orthogonal, unitary and symplectic 
symmetry, characterizing the energy levels of an ergodic 
particle without magnetic field, with magnetic field, and 
with spin- orbit- scattering, respectively ]l7| |. The con- 
fusing fact, that random orthogonal matrices result for 
the ASD in a functional integral over compact symplectic 
Q- matrices and vice versa, results from the sign change 
of the Grassman variables under time reversal. 

We can derive the corresponding Hamiltonian H by 
means of the transfer matrix method, reducing the one- 
dimensional integral in Eq. (3) to a single functional 
integral. Thus, the ASD is obtained in the simple form 
of Eq. (H|), with the effective Hamiltonian 

H{u = 0) = -J_(-4A g - ^X 2 Tr Q [Q,r,] 2 ). (5) 

Aq is that part of the Laplacian on the symmetric space, 

which couples to Tr[A 3 Q]. X = 2naL C u(y 2 ~) il/2) B/<po, 
where 4>q = q/h is the flux quantum. 

The problem is now equivalent to a particle with 
"mass" (a/8)Lcu{E) moving on the symmetric space of 
Q in a harmonic potential with "frequency" X/(aLcu)i 
and in an external field «a(7r/4)w/(LA), in "time" x, the 
coordinate along the wire. To find the ASD as a function 
of uj and the length of the wire L, one can do a Fourier 
analysis in terms of the spectrum and eigenfunctions of 
the effective Hamiltonian at zero frequency, H(uj = 0) 
Pq| . For L 3> £ c , only the gap between the ground and 
the first excited state of H(cu = 0), determines the ASD. 

Without magnetic field, B — 0, the Laplacian is 



A Q =d Xc (l~\ 2 c )dx c +2 



l-X 



Ac 



Ac: 



2 



4-^(1 - X 2 D )9x D , 



(6) 



where Xc.d € [ — 1, 1] . Its ground state is 1 and the first 
excited state is Ac Ad- Thus, the gap is 



E G (B = 0) = 16/Lcu- 



(7) 



For moderate magnetic field, with the condition 
Lcuiy 2 )'" 1 ^ 2 ^ B 3> 4>q — h/q, all degrees of freedom aris- 
ing from time reversal invariance are frozen out, due to 
the term Ttq [Q, t 3 ] 2 = 16(A C - 1) which fixes A c = 1. 
Then, the Laplacian reduces to 



A Q =d XD (l->?D)dx D . 



(8) 



Its eigenf unctions are the Legendre polynomials. There 
is a gap above the isotropic ground state of magnitude 



E G {X > 1) = 8/L 



cu- 



(9) 



For moderate magnetic impurity scattering exceeding the 
local level spacing, 1/t s > Ac, the Laplacian is given 
by Eq.(8). Thus, due to a = 2, the gap is reduced to 
Eg(1/ts > Ac) = i/Lcu- For moderately strong spin- 
orbit scattering 1/tso > Ac, the Laplace operator is 



A, 



E 

1=1,2 



d Xl (l-X 2 )d x 



(10) 



where Ai,2 6 [—1, 1]. The ground state is ijjQ = L the first 
excited state is doubly degenerate, tpu — Ai, V12 = A 2 . 
Thus, the gap is the same as for magnetic impurities, 



E g {1/tso > Ac) = 4/L 



cu- 



(11) 



An external magnetic field lifts this degeneracy but does 
not change the gap. 

Using the crossover in energy level statistics as the def- 
inition of a localization length as above, we get in a quasi- 
1 -dim. wire, 



£ c = 1/E G (0) = (l/W)PL cu , 



(12) 



where (3 = 1,2,4 corresponding to no magnetic field, fi- 
nite magnetic field, and strong spin- orbit scattering or 
magnetic impurities, respectively. Comparing with the 
known equation for the localization length, L Cl we find 
that the dependence of the ratios (3 on the symmetry are 
in perfect agreement with the result as obtained from the 
spatial decay of the density- density- correlation function 
0, while it defers by the overall constant 1/8. 

This relation can be proven directly. The ASD at zero 
frequency C(0)l of the wire of length L, becomes, when 
the wire is divided into two parts, C(Q) 2 L , 2 . For L — > 00, 
we find that the relative difference is: 



f(L) = 



C(0) 2 



L/2 



C(0)i 



-l = 2exp(- J LS G /2), 



(13) 



exponentially decaying with the length L. f(L) can be 
estimated, following an argument by Mott M): When 
the two halves of the wire get connected, the Eigen- 
states of the two separate halves become hybridized 
and the Eigenenergy of a state ip n is changed by 
±Ac exp(— 2x n /Lc)- x n is random, depending on the 
position of an eigenstate with closest energy in the other 
half of the wire. Thus, averaging over x n gives: 



f(L) ~ + exp(-4L/L, 



c 



(14) 



Comparison with Eq. ( |l3| ) yields indeed 1/Lc = 8E G . 
The Crossover Behaviour of the Localization Length 
This direct relation of the ASD to the spatial overlap 
between wavefunctions allows us to extend this approach 
to get an analytical solution for the crossover behaviour 
of the localization length and the local level spacing as 
a magnetic field is turned on, and without strong spin- 
orbit scattering. While a heuristic approach g and nu- 
merical studies J9| seemed to indicate a continous increase 
of the localization length, the analytical result [jl0| does 
show that both limiting localization lengths L c (f3 = 1) 
and L c ((3 — 2) are present in the crossover regime and 
that there is no single parameter scaling. This is ex- 
plained by arguing that the far tails of the wavefunc- 
tions do cover a large enough area to have fully broken 
time reversal symmetry, decaying with the length scale 
L C (B = 2) even if the magnetic field is too weak to af- 
fect the properties of the bulk of the wavefunction, which 
does decay at smaller length scales with the shorter lo- 
calization length L c {j3 — 1), corresponding to the time 
reversal symmetric case. The quantity studied there is 
the imprurity averaged correlation function of local wave- 
function amplitudes and its momenta at a fixed energy 
e: Y(e) =< J2 a | ^a(O) | 2 | | 2 S(e - e a ) >. It is 

averaged over a distribution of eigenfunctions in different 
impurity representations. Thus, each eigenfunction could 
decay exponentially with a single localization length, but 
with a distribution which has two maxima, at L c (f3 — 1) 
and L c {(3 — 2), whose weight is a function of the mag- 
netic field in the crossover regime. It is interesting to ask 
if this property of the eigenfunctions does have a conse- 
quence on the energy level statistics as well, or, if that is 
still governed by a single parameter as the magnetic field 
is varied. 

The effective Hamiltonian for moderate magnetic fields 
is given, without spin dependent scattering, a = 1, by: 



H 



Leu 



(-4A Q + X 2 (l-\ 2 C )), 



(15) 



where the Laplacian is Eq.|| and X = 2ir(j>/4>o, where <j> 
is the magnetic flux through an area A — Lct/(y 2 )^l/2). 

In the limit X — > the ground state and first excited 
state approach 1,AcAd, respectively. In the limit X 3> 
1, X c becomes fixed to 1. Thus, the Ansatz -00 (Ac) = 
exp(A)(l - Xc))' and ^i(Ac,Ad) = A c Ad exp(Ai(l - 



3 



Xq)), solves Hip 



Eip to first order in z = X 2 (l 



A 2 C ) 



and yields the gap: 

E G {X) = 4(2 + ^49 + X 2 - ^25 + X 2 )/L cu . (16) 

This solution is valid in both the limits X <C 1 and 
X ^> 1, interpolating the region X « 1. The result- 
ing ratio of local energy level spacings Ac{B)/Ac(0) = 
Eg(B)/ Ec(0), is shown in Fig. 5 to be in excellent quali- 
tative agreement with experimental data for the magnetic 
field dependent activation energy measured recently in 
transport experiments [Q. 

There is a quantitative discrepancy by a factor .6 be- 
tween the best fit X = .036H/Oe, and X = 27r0/</>o, 
4> = fioH Lcu(y 2 ) < " 1 ^ 2 ' > ■ With the experimental param- 
eters a — I, Leo — .61/zm, width W = .2/im of sam- 
ple 5 in Ref. and y 2 = W 2 /12 for a 2- dimensional 
wire, it yields X = .021H/Oe. We note that smooth 
confinement can give y 2 > W 2 /12, thus explaining this 
discrepancy and the observed difference between W as 
obtained from the sample resistance and estimated from 
the analysis of the weak localization magnetoresistance, 
which also depends on y 2 plj ). When W 2 B > 0o or 
H > H* = 2.1 * 10- n Oe[mVtV 2 ] = 5250e, the deriva- 
tion has to be extended to the 2- dimensional nonlinear 
sigma- model, and Eq. [r] seizes to be valid. Then, Lc 
increases as function of magnetic field more strongly. 

The asymptotic behaviour, of 5Lc(B) ~ B 2 for small 
and ~ 1/B at large magnetic fields of Eq. (|l^) does agree 
with the suggestion by Bouchaud || . 
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FIG. 1. The activation gap ratio To(H) /Tb(0) as a func- 
tion of the magnetic field H in Oe of sample 5 measured at 
temperature T = .3K as reported in Ref. 1, together with 
the theoretical curve. 



The two length scale physics of Ref. |T(| has thus no 
consequence for the energy level statistics as studied with 
the ASD. Since we could also show a direct relation be- 
tween this spectral statistics and the exponential decay 
of localized eigenf unctions, it is suggestive to be the non 



self averaging property of the ASD which washes out the 
two scale physics. 
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